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ABSTRACT 

We have undertaken deep optical imaging observations of three 6.2 < z < 6.5 quasar fields in the i' and z' 
filters. These data are used to search for foreground galaxies which are gravitationally lensing the quasars and 
distant galaxies physically associated with the quasars. Eoreground galaxies are found closer than 5" from the 
lines-of-sight of two of the three quasars. However, the faintness of these galaxies suggests they have fairly 
low masses and provide only weak magnifications (/i ^ 1.1). No convincing galaxies physically associated with 
the quasars are found and the number of /'-band dropouts is consistent with that found in random fields. We 
consider the expected dark matter halo masses which host these quasars under the assumption that a correlation 
between black hole mass and dark matter halo mass exists. We show that the steepness of the high-mass tail 
of the halo mass function at this redshift, combined with realistic amounts of scatter in this correlation, lead 
to expected halo masses substantially lower than previously believed. This analysis can explain the lack of 
companion galaxies found here and the low dynamical mass recently published for one of the quasars. 

Subject headings: cosmology: observation - gravitational lensing - quasars: general 


1. INTRODUCTION 

Active supermassive black holes provide a useful means of 
locating very distant, massive galaxies. Studies of the host 
galaxies of luminous active galactic nuclei (AGN) at redshifts 
up to z w 2 show they are associated with galaxies with lumi¬ 
nous rest-frame-optical stellar populations, corresponding to 
(Kukula et al. 2001; Ridgway et al. 2001; Hutchings 
et al. 2002; Willott et al. 2003a; Dunlop et al. 2003). Eur- 
thermore, the strong coiTelation between black hole mass and 
stellar bulge luminosity/mass observed locally (Magorrian et 
al. 1998; Gebhardt et al. 2000; Eerrarese & Merritt 2000) 
shows that the most massive galaxies are those whose black 
holes have accreted the greatest mass and therefore were the 
most luminous AGN. 

It is now possible to discover luminous quasars out to a 
redshift of z = 6.4. Sources at this redshift are observed as 
they were w 13 Gyr ago (93% of the age of the universe). 
The Sloan Digital Sky Survey (SDSS) has made a spectac¬ 
ular breakthrough in locating such quasars and now 12 are 
known at z > 5.7 (Ean et al. 2004). Black hole masses in these 
quasars are usually estimated by making the assumption that 
the quasars are accreting at the Eddington limit (e.g. Ean et 
al. 2001). This method is consistent with a black hole mass 
measurement from the kinematics of the broad emission line 
gas (Willott, McLure & Jarvis 2003b). Under the Eddington 
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argument, the black hole masses of these quasars are mostly 
in the range 1 - 5 x IO^Mq. These black hole masses could 
be much lower in the case of strong gravitational lensing or 
beaming. Observations so far have not shown such effects 
(Ean et al. 2003; Richards et al. 2004; Willott et al. 2003b). 

The space density of these luminous SDSS quasars is ex¬ 
tremely low, p = (6.4 ±2.4) X 10“'°Mpc“^ (Ean et al. 2004). 
The rarity of these objects and their large black hole masses 
makes it tempting to associate these quasars with the rarest 
peaks in the dark matter density distribution (Ean et al. 2001). 
The dark matter halos hosting the quasars would, under this 
hypothesis, have a mass of > IO'^Mq. Such halos would 
merge with smaller halos and in the present day would be 
identified as massive galaxy clusters with a giant elliptical 
galaxy hosting the dormant supermassive black hole at the 
bottom of the potential well. However, it is also possible that 
the SDSS quasars reside in much more common, lower-mass 
halos and that the correlation between halo mass and black 
hole mass is not well established at this early epoch. Deter¬ 
mining the masses of the quasar host halos is extremely im¬ 
portant for understanding the connection between black hole 
and galaxy growth in the early universe. 

The evidence from sub-millimetre observations paints a 
rather confusing picture. On the one hand, thermal emis¬ 
sion from dust has been detected from several of the quasars. 
These detections imply huge dust masses and star formation 
rates of several thousand Mq yr“* (Bertoldi et al. 2003a; Prid- 
dey et al. 2003). At this rate, the stellar mass of an L* ellipti¬ 
cal could be built up in ~ 0.1 Gyr. The clustering of galaxies 
with such high star formation rates at lower redshift (z ~ 3) 
is not well constrained, but early results suggest these galax¬ 
ies are associated with massive dark matter halos (Blain et al. 
2004, although see Adelberger 2005). One of these quasars, 
SDSS J11483-5251, has also been detected in molecular car¬ 
bon monoxide transitions (Bertoldi et al. 2003b; Walter et 
al. 2003). The CO spectra have a relatively narrow velocity 
profile with width 280 kms“'. Walter et al. (2004) obtained 
high-resolution imaging of the CO emission and confirm it 
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comes from a compact structure (a few kpc). These obser¬ 
vations allow a dynamical mass estimate for the mass within 
the central 2.5 kpc of the galaxy. This mass is comparable to 
the inferred molecular gas mass and is an order of magnitude 
lower than the mass predicted by assuming that this quasar 
resides in one of the rare peaks corresponding to dark matter 
halos of > IO^^Mq. Although there are some uncertainties in 
the dynamical mass estimate, particularly the geometry and 
inclination of the gas, this result casts serious doubt on the be¬ 
lief that the SDSS quasars pinpoint the most massive galaxies 
at high redshift. 

Another approach to determining the mass of the quasar 
host dark matter halo is via a search for companion galaxies. 
If the SDSS quasars are formed in the rarest density peaks, 
then these correspond to large-scale overdense regions and 
hence the number of dark matter halos located nearby is sub¬ 
stantially above the cosmic mean (Kaiser 1984; Barkana & 
Loeb 2004). Hence the rarest peaks at z Ri 6 are likely to be 
the sites of the first proto-clusters. If the quasars really do oc¬ 
cupy very massive halos, we would expect to see star forming 
galaxies in their vicinity. To attempt to find these compan¬ 
ion galaxies we have carried out deep optical imaging around 
the three highest redshift quasars from the sample of Fan et 
al. (2003). In this paper we present these data, describe a 
search for foreground galaxies which may be gravitationally 
lensing the quasars and a search for star forming galaxies at 
the quasar redshift. Finally, we discuss the implications for 
the host halos of the most distant quasars. Cosmological pa¬ 
rameters of Ho = 70kms“* Mpc“\ Hm = 0.3 and Ha = 0.7 are 
assumed throughout. 

2. OBSERVATIONS 

We present observations for three quasars: SDSS 
J103027.10H-052455.0 (SDSS J1030 h-0524; z = 6.30), SDSS 
J104845.05H-463718.3 (SDSS J1048h-4637; z = 6.20) and 
SDSS J114816.64H-525150.3 (SDSS J1 148 h-5251; z = 6.42). 
The three quasar fields were imaged using the GMOS- 
North imaging spectrograph on the Gemini-North Telescope. 
GMOS-North uses an array of three 2048 x 4608 pixel EEV 
CCD detectors. The pixel scale is 0.073" per pixel giving a 
useful field-of-view for imaging of 5.5' x 5.5'. Since the typ¬ 
ical seeing size of the observations is in the range 0.5 to 0.7", 
the pixels were binned by a factor of 2 in both directions to 
speed up readout and data processing. The observations were 
carried out in queue mode during November and December 
2003. Typical exposure times are « 2 hours in the z'-band 
and Rs 3 hours in the /'-band. The relative exposure times 
were designed to give similar sensitivity in the two bands for 
very red objects with colours of /'-z' ~ 1.5. More details of 
the observations are given in Tabled 

The GMOS-North detectors have significant fringing at the 
red end of the optical wavelength regime. This is of particu¬ 
lar concern for the observations presented here which use the 
red /' and z' filters. The fringe pattern remains constant over 
time, so the fringes can be successfully removed if one ob¬ 
tains sufficient data to construct a high signal-to-noise fringe 
frame free of stars and galaxies. To enable the construction of 
these fringe frames, the observations of each field in each fil¬ 
ter were split into 20 to 30 individual frames and the telescope 
was dithered in a 3 x 3 pattern with offsets of 10" in each di¬ 
rection. Median-combining all the data in each filter from the 
three fields enabled the construction of fringe frames free of 
astronomical objects. 

Reduction of the imaging data was carried out using stan- 


TABLE 1. GMOS-N IMAGING OBSERVATIONS. 


Quasar 

Band 

Exposure 
time (s) 

3cr limiting 
mag (AB) 

Seeing 

(") 

SDSS J1030+0524 

z' 

8100 

26.3 

0.68 

SDSSJ1048+4637 

z! 

5850 

26.2 

0.61 

SDSS J1148+5251 

z! 

7650 

26.2 

0.65 

SDSSJ1030+0524 

i' 

9900 

27.5 

0.65 

SDSSJ1048+4637 

i' 

12150 

27.7 

0.54 

SDSS 11148+5251 

i' 

13050 

27.6 

0.67 


dard procedures. Most of the reductions were performed us¬ 
ing tasks in the IRAE Gemini package which were specif¬ 
ically designed for GMOS. The first step is removal of the 
bias level using a bias frame constructed from many (> 10) 
bias observations (with the same binning) carried out during 
the same month. The images were flat-fielded using flat-held 
frames generated from observations of the twilight sky. The 
three separate CCD images were then mosaiced together into 
one large image. The fringe frames generated from all ob¬ 
servations within each filter were then scaled and subtracted 
from each image. All images were then inspected and in some 
cases a different scaling was applied to better subtract off the 
fringes. The images were then scaled to correct for atmo¬ 
spheric extinction. 

All the images in each Alter of the same field were then po¬ 
sitionally registered using detected objects and combined into 
one image rejecting pixels deviating by more than 2.5 sigma 
from the median and using a bad pixel mask. This method 
successfully removed CCD defects and cosmic rays whilst 
not rejecting counts from actual objects. The combined im¬ 
ages still showed a low-level, large-scale, residual background 
which varied across the images. The large-scale background 
was successfully removed by applying the background sub¬ 
traction method of the Sextractor software (Bertin & Arnouts 
1996) with a mesh size of 19". The /' images for each field 
were then shifted so that objects would appear in the same 
locations in the /' and z' images. Astrometry was performed 
by using the known geometric distortion of GMOS-North on 
the sky and setting the positions of the quasar targets to the 
positions given in the SDSS survey. Photometric calibration 
was achieved using photometric standard stars observed dur¬ 
ing the same nights as some of our observations. 

3. OBJECT DETECTION, PHOTOMETRY AND COMPLETENESS 

Detection of objects in the images was performed using the 
Sextractor software. The z'-band was selected as the primary 
detection waveband since z > 6 galaxies are expected to have 
/'-z' > 2 and may therefore be undetected at /'. Sextractor 
was run in “double-image” mode to determine /'-band mea¬ 
surements for objects detected in z'. The edges of the im¬ 
ages do not contain data at all dither positions and hence have 
lower sensitivity, contain artifacts and have a varying back¬ 
ground. Objects in these regions were excluded from the ob¬ 
ject catalogues. The useful area of the catalogues is 27, 27, 28 
square arcmin for SDSS J1030H-0524, SDSS J1048H-4637 and 
SDSS J1148H-525, respectively. 

Magnitudes on the AB system were measured in circular 
apertures of diameter 1.5". This size was chosen because it is 
greater than twice the seeing size, is much greater than the size 
scales of known z > 6 galaxies (Bouwens et al. 2004) and has 
lower magnitude errors than larger apertures. Aperture cor¬ 
rections were applied statistically to the z'-band magnitudes 
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Fig. 1.— Binned number counts in the z'-band in the fields of the three 
z > 6 quasars. The z^-band magnitudes have been aperture-corrected using 
the prescription given in Sec.|^ Error bars include only Poisson errors. The 
solid curve is the mean counts averaged over the three fields. The open stars 
show the z^-band number counts from the deeper and wider-area Subaru Deep 
Field (Kashikawa et al. 2004). The counts are consistent given the size of the 
error bars up to t! =26. Beyond ^ — 26 the quasar field counts turn over 
indicating severe incompleteness. 


by fitting a linear function to the difference between the total 
magnitude and aperture magnitude as a function of aperture 
magnitude. The best-fit relation is Ztotai“^ap = 3.58-0.133z'p 
which gives an aperture correction of 0.25 mag at tI = 25. 
Note that this relation will lead to large aperture corrections at 
bright magnitudes, which may be inappropriate for compact 
sources such as stars, but we are in general only interested in 
faint objects for which the aperture corrections are reasonable 
(z' > 23). A similar procedure was not adopted for the Z'- 
band magnitudes since the uncorrected aperture magnitudes 
are used to measure the i' -z' colours. 

The rms noise in the sky background was measured to de¬ 
termine the magnitude limits of the images. Magnitude limits 
are quoted as 3 ct limits in 1.5" apertures. Typical magnitude 
limits are z! ~ 26.2 and /' k, 27.6 (see Table[0. The relative 
depths of the images are suitable for detecting very red objects 
with i' -zl >1.5. 

To assess the completeness of the z'-band catalogues we 
consider both the observed number counts and the recovery 
of simulated objects. Binned number counts for the z'-band 
images for all three fields are shown in Fig. [2 The counts 
in the three fields do not differ significantly from each other. 
They agree well with the z'-band counts determined from a 
deeper and much larger area survey (0.2 square degrees or 
thirty times the GMOS-North field-of-view) of the Subaru 
Deep Field by Kashikawa et al. (2004). The number counts 
in the quasar fields begin to change slope at z' > 25.5 and turn 
over at z! = 26, indicating this is where the sample becomes 
incomplete. 

The source recovery as a function of magnitude was de¬ 
termined by populating the images with artificial galaxies 
and then using Sextractor to attempt to detect these objects. 



z' mag (AB) 

Fig. 2.— Completeness ratio against aperture-corrected z'-band magnitude 
derived from recovery of simulated galaxies as detailed in Sec.|3 The curves 
are quite similar for the three different fields. Dotted lines indicate the loca¬ 
tion of a completeness ratio of 0.8 and the adopted complete magnitude limit 
of z' = 25.5. 


About 10000 artificial galaxies with magnitudes in the range 
23 < z' < 27 were placed into copies of the z' images of each 
quasar. Regions of the images already occupied by objects 
were masked out of the process to eliminate incompleteness 
due to blending. Sextractor was run twice; firstly on images 
containing only the artificial galaxies and a very low noise 
level and then on the actual quasar field images with the ar¬ 
tificial galaxies inserted. The ratio of the number of artificial 
objects detected in the quasar field images to the number in 
the low noise images gives the completeness. This complete¬ 
ness ratio is plotted as a function of magnitude for the three 
quasar fields in Fig.|2l The completeness in all the fields is 
fairly flat at close to 1 up to z' = 25.2 and then begins to de¬ 
cline. The rapid decline occurs at z' > 25.5 and the complete¬ 
ness drops to 0.5 by z = 26.0. All the fields have completeness 
> 0.8 at z' = 25.5 and we adopt this as the magnitude at which 
completeness begins to become an issue. This analysis with 
simulated objects agrees well with the results for the number 
counts discussed previously. 

4. SEARCH FOR FOREGROUND GALAXIES TENSING THE 
QUASARS 

The combined effects of a steep luminosity function and 
the high optical depth to z ~ 6 mean that gravitational lens- 
ing is expected to be particularly important for surveys of lu¬ 
minous quasars at high redshift, such as the SDSS z > 5.7 
quasar survey. Predictions for the fraction of multiply im¬ 
aged quasars in the SDSS for various forms of the luminosity 
function have been made by Wyithe & Loeb (2002a;b) and 
Comerford, Haiman & Schaye (2002). Fan et al. (2003) and 
Richards et al. (2004) discuss high-resolution imaging (rang¬ 
ing from 0.1" to 0.8") of all seven known z > 5.7 quasars at 
the time and found that none of them appears to be multiply- 
imaged. However, the lack of multiple images of the SDSS 
quasars does not necessarily mean that lensing is unimpor- 
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Fig. 3.— CMOS images in the i'- and z'-bands centred on the three 
quasars SDSS J1030-I-0524, SDSS J1148-1-5251 and SDSS J1048-H4637 (left 
to right). The final column shows the images of SDSS J1048-I-4637 after sub¬ 
tracting off a scaled point-spread function at the location of the quasar to 
show more clearly the nearby galaxies. There are some residuals from the 
psf-subtraction process which tests on stars in the field show are likely to be 
artifacts. Each box is 10" on a side, centred on the quasar with north up and 
east to the left. Contour plots are shown with a logarithmic scale at approx¬ 
imately 1.5,3,5,9,17,31,56 x the background rms per pixel. Greyscale plots 
have a linear stretch ranging from -0.5 X the background rms (white) to +4 X 
the background rms (black). 


tant. For simple singular isothermal sphere models the mag¬ 
nification is p, < 2 if the source is singly-imaged. More re¬ 
alistic potentials including isothermal ellipsoids, cluster-scale 
halos and micro-lensing can give rise to higher magnifications 
for singly-imaged quasars (Keeton, Kuhlen & Haiman 2004; 
Wyithe & Loeb 2002a). 

Shioya et al. (2002) found a faint galaxy close to the line- 
of-sight of the z = 5.74 quasar SDSS J1044-0125. The lens- 
ing galaxy has a magnitude i'(AB) = 24.3, separation from the 
quasar 0=1.9” and a likely redshift in the range 1.5 < z < 2.5. 
Spurred on by this discovery, this group observed two more 
high-redshift quasar fields (SDSS J1030-I-0524 at z = 6.30 and 
SDSS J1306-1-0356 at z = 5.99), but failed to find galaxies pro¬ 
jected closer than 6" and 3" to these quasars, respectively (Ya- 
mada et al. 2003). 

We have imaged three quasar fields and in this section we 
use our data to search for foreground galaxies which may be 
lensing the quasars. One of these fields, SDSS J1030H-0524, 
has been imaged before by Yamada et al. (2003), but our 
images are 1-2 mags deeper and have better resolution than 
those in Yamada et al. and Shioya et al. Images of the regions 
surrounding the quasars are shown in Fig.|3 Contours and 
greyscales are shown for the same fields to highlight different 
aspects of these high dynamic range images. The contours 
show smooth isophotes in the quasar flux which rule out the 
existence of very close, bright companion galaxies and grav¬ 
itational lensing image splitting in the quasars at the limit of 
our resolution. The ellipticities and FWHM of the quasars are 


Table 2. Objects detected within 5" of the quasars. 


Object 

Object position (J2000) 

0/" 

/' (AB) 

z! (AB) 

a 

10:30:27.26 -1-05:24:55.3 

2.5 

27.54 ±0.35 

>26.3 

b 

10:30:27.08 -1-05:24:58.2 

3.5 

26.69 ±0.17 

26.12 ±0.28 

c 

10:30:26.99 -1-05:24:50.6 

4.6 

26.05 ±0.10 

>26.3 

d 

10:48:45.19 -1-46:37:17.4 

1.7 

25.67 ±0.05 

24.97 ±0.11 

e 

10:48:45.10-1-46:37:16.1 

2.2 

27.32 ±0.22 

> 26.2 

f 

10:48:44.92 -1-46:37:15.9 

2.7 

25.37 ±0.04 

24.82 ±0.10 


Note. — Objects detected at /'-band within a 5" radius of the quasars. 
Three objects were detected in each of the fields of SDSS J1030-I-0524 and 
SDSSJ1048-I-4637. Nothing was detected wi t hin 5" of SDSS J1148-H5251. 
Object refers to the labels in Fig.|^ 9 is the separation of the object and the 
quasar. Magnitudes are given as observed in a 1.5" diameter aperture. Limits 
at z'-band are 3(7 in the same size aperture. 


consistent with those of stars in the images. 

We searched for objects close to the quasars by inspection 
of smoothed images. We only consider objects within a 5" 
radius of the quasars, since lensing by galaxy-sized masses 
is only effective with impact parameters of a few arcsec. In 
addition, the probability of finding a galaxy at random within 
a circle of 5" radius is large at the magnitude limit of the i'- 
band images (this probability is 0.9 for i' < 27). 

Three galaxies were detected in the fields of each of 
SDSSJ1030H-0524 (labelled a, b, c) and SDSS J1048H-4637 
(d, e, f) and none in the field of SDSS J11483-5251. The 
details of these galaxies are given in Table|2] All of these 
galaxies are visible in the /'-band and their /' - z' colours (or 
colour limits) show they are foreground to the quasars. For 
SDSS J1030H-0524, the objects are all rather faint (/' > 26) 
and distant from the quasar (> 2.5"). For SDSS J1048H-4637, 
the closest galaxy (d) is only 1.7" from the quasar. To mea¬ 
sure accurate photometry for this galaxy it was necessary to 
subtract the bright unresolved quasar from the image (far-right 
panels of Fig.|3l. Galaxy d has 1.5" aperture magnitudes of 
/' = 25.67 ± 0.05, z' = 24.97 ± 0.11. The galaxy 2.7" from the 
quasar (f) has similar magnitudes. Both these galaxies have 
fairly red colours with i' -z' = 0.70 and 0.55, respectively. 
Only 20% of galaxies in our images with similar magnitudes 
have i'-z' > 0.5. 

4.1. Lensing galaxy masses, redshifts and magnifications 

In order to determine the lensing magnification due to such 
galaxies one needs to know the galaxy mass (or equivalently 
the velocity dispersion) and redshift. Under the simplify¬ 
ing assumption that the lens potential can be modelled as an 
isothermal sphere the magnification /r is p, = 0/{6-9^), where 
0 is the angle from the lens galaxy to the source, is the Ein¬ 
stein radius given by 0^ = 47r(cry/c)^(DLs/Dos), where Dls 
is the lens-source angular diameter distance and Dos is the 
observer-source angular diameter distance. 

Neither the masses nor redshifts are strongly constrained 
by our photometric observations, but we will use the infor¬ 
mation available to determine likely values. The red colour 
of i' —z' = 0.70 for galaxy d close to SDSS J1048H-4637 is not 
consistent with the colour of known galaxies at low redshifts. 
Therefore we can at least put a lower limit on the redshift 
of this galaxy. Also, galaxy colours at high redshifts are ex¬ 
pected to be blue in the absence of dust reddening, due to 
the prominence of young stellar populations. We used the 
HyperZ code (Bolzonella, Miralles & Pello 2000) with the 
2 photometric points to determine if this colour leads to any 
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constraints on the redshift. Using a range of evolving galaxy 
templates with moderate dust reddening (0 < Ay < 0.5) we 
find a 90% confidence range for the redshift is 0.4 < z < 2.1. 
There is also a secondary solution at 5.1 < z < 5.5 where the 
red colour is due to the Lyman break entering the i'-band. 
Considering the redshift distribution of galaxies of this magni¬ 
tude it is most plausible that the galaxy is at the lower redshift 
range. Note that if this galaxy is highly reddened by dust, then 
the redshift is no longer well-constrained. 

The mass or stellar velocity dispersion of the galaxy can 
be estimated from the photometry by assuming a redshift and 
star formation history. Given the red colour of this galaxy 
we assume that it is an early-type which has undergone pas¬ 
sive evolution after a starburst at high redshift and use the 
evolving elliptical galaxy template from Bruzual & Chariot 
(2003) which formed at z = 7 (the exact formation redshift as¬ 
sumed does not critically change the results presented here). 
The correlation between i'-band absolute magnitude and 
observed at low redshifts (Bernard! et al. 2003) was used to 
relate the evolved absolute magnitude to cr^. For three possi¬ 
ble lens redshifts of z = 0.4,1,2 the calculated stellar velocity 
dispersion is dv = 55,87,127kms“'. Note that the Bernard! et 
al. relation is only derived for z < 0.3 and cr,, > lOOkms”', 
so its use here involves extrapolation in both redshift and ve¬ 
locity dispersion. Given the redshift of the source is z = 6.20 
and 9 = 1.7" this leads to Einstein radii 9e = 0.07,0.13,0.18" 
and magnifications p = 1.04,1.08,1.12 for lens redshifts of 
z = 0.4,1,2, respectively. 

It is clear that, whatever the redshift of this lensing galaxy, it 
is not very massive and provides only a small magnification to 
the flux of SDSS J1048H-4637. Therefore it makes a negligi¬ 
ble difference to the derived luminosity and black hole mass. 
Given that the closest galaxy to SDSS J1030H-0524 (galaxy a) 
is even fainter and further away, the magnification in that case 
will be even lower. An increase in the magnification could 
come about if (i) the galaxy lies in a moderately rich high- 
redshift cluster; or (ii) the galaxy resides in a dark matter halo 
with higher than average mass-to-light ratio. 

For SDSS J1044-0125, Shioya et al. (2002) estimated the 
lens galaxy to have a velocity dispersion of cr,, ~ 200 km s“' 
and be located at redshift z « 2. According to Shioya et al. 
this would give a magnification of p = 2. However we have re¬ 
peated their calculations and find their values for the Einstein 
radius and magnification to be erroneous (this error is also 
present in Yamada et al. 2003). For z = 2 and = 200km s“\ 
the magnification is actually fjL = 1.3. The range of plausi¬ 
ble lens redshift and velocity dispersion combinations given 
in their paper lead to a range ofl.l<p<1.5. As with 
SDSS J1048H-4637, the corrections to the derived luminosity 
and black hole mass of SDSS J1044-0125 are small. 

5. SEARCH FOR Z > 6 GALAXIES IN THE QUASAR EIELDS 

Galaxies at redshifts z > 5.7 can be identified by the sharp 
drop in flux across the Lyman break leading to a very red 
i' - z' colour. The z > 6 SDSS quasars have z' ~ 20 and 
colours measured from our images of i' -z' = 3.25,2.98,3.25 
respectively for SDSS J1030H-0524, SDSS J1048H-4637 and 
SDSS J1148H-525L The spectra of quasars and star form¬ 
ing galaxies over the rest-frame wavelength range 90-140 nm 
probed by the i' and z' filters are dominated by a large break 
due to absorption by neutral hydrogen. Therefore one would 
expect companion galaxies to have comparable i' — z' colours 
to the quasars. The fact that the light from a companion 
galaxy would pass through IGM neighbouring that which 


z mag (AB) 

21 22 23 24 25 26 27 



Eig. 4.— Colour-magnitude (aperture-corrected) diagrams for objects de¬ 
tected in z' in the quasar fields. The completeness limit of z' = 25.5 is 
shown with a dot-dashed line. The high-redshift galaxy selection criterion 
of i' — z' > 1.5 is marked with the dashed line. The dotted line represents 
the colour of an object which is just detected at the 3cr level in the f-band 
as a function of the z'-band magnitude. Most objects with measured zero or 
negative flux in the /'-band are plotted at /'—z' = 3. The exceptions to this are 
objects which pass the /'-band dropout selection criteria discussed in Sec.|^ 
These are shown with star symbols and lower limits on the /'-band 3 it line. 
The dots at /' -z' >1.5 indicate sources which are detected at less than Act 
at z'-band and hence have very large uncertainties on their magnitudes and 
colours and in some cases may be spurious. The labelled curves show the 
colour and magnitude as a function of redshift for an evolving L* elliptical 
(upper curve) and a non-evolving L* Lyman-break galaxy (lower curve) - see 
text for more details. 


the quasar light passes through, further strengthens the idea 
that the colours of companion galaxies are expected to be 
i'-z'-3. 

Simulated colour-magnitude tracks for two different types 
of galaxy as a function of redshift are shown on Fig.0 Model 
galaxy spectra were generated from the Bruzual & Chariot 
(2003) spectral synthesis code with Lyman forest absorption 
evolution matching the observations of Songaila & Cowie 
(2002). The upper curve is a present-day L* elliptical which 
formed all of its stars in a starburst starting at z = 10 with a 
characteristic timescale of 1 Gyr and evolved since without 
merging. There is no dust extinction assumed for this model, 
but in reality the dust extinction would increase with redshift 
(due to evolution and k-correction) making the galaxy redder 
and fainter than plotted at higher redshifts. The lower curve 
is a L* Lyman-break galaxy model where the galaxy is ob¬ 
served 0.5 Gyr into a constant star formation rate starburst. It 
is clear that the only possible low-redshift contaminants with 
i' -z' >1.5 are galaxies with a high level of dust reddening. 
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Also plotted on Fig.|4]are colour-magnitude diagrams con¬ 
structed from the z'-band-selected catalogues described in 
Sec.|3for each of the three quasar helds. Most objects have 
colours in the range 0 < i' -z' < 1 as is well known from pre¬ 
vious surveys (Dickinson et al. 2004; Capak et al. 2004). 
These colours are in agreement with the model galaxy tracks 
at z < 5. 

A search was made for objects which could plausibly be 
high-redshift galaxies. The /'-band dropout selection criteria 
adopted were signal-to-noise ratio in the z'-band > 4 and a 
colour of i'-z' > 1.5. Possible candidates were inspected and 
magnitudes checked to ensure their unusual colours are not 
spurious. A total of three objects satisfying these criteria were 
found - these are shown with hlled symbols on Fig.|3 These 
all have magnitudes in the range 25 < z' < 26 and only one 
of them is detected at the 3cr level at /'-band. This object has 
z' = 25.10 and a colour of /'-z' = 1.56 ±0.22. The other two 
objects have z' « 25.6 and only lower limits on their /'-z' 
colours which lie in the range 1.7 - 1.9. 

In Sec.|3we showed that our z'-band images are complete 
at the greater than 80% level to z' = 25.5. Fig.|4]shows that at 
z' < 25.5, every single z'-band object detected on our images 
has a counterpart at the > 3(7 level in the /'-band image. At 
fainter magnitudes, this is no longer true and our constraints 
on the /' - z' colours of objects at these magnitudes becomes 
very weak due to uncertainty in both the z' and /' magnitudes. 
Therefore we will limit further analysis to objects brighter 
than the magnitude limit of z' = 25.5. 

At z' < 25.5 we hnd one /'-band dropout in the three quasar 
helds. This object has /'-z' = 1.56 ±0.22 which is quite close 
to the colour selection value and the size of the uncertainty 
means it is quite plausible that photometric errors have scat¬ 
tered the colour into the dropout range. The measured colour 
is about 7(7 away from the colours of the SDSS quasars, which 
suggests that if it is a high-redshift galaxy, then it is most 
likely foreground to the quasars with a redshift in the range 
5.7 < z < 6. The galaxy is located a projected distance of 
81" from the quasar (approximately halfway from the quasar 
to the edge of the GMOS field). Therefore we conclude that 
there are no plausible galaxies brighter than z' = 25.5 associ¬ 
ated with any of these three quasars. 

We now consider the number of /'-band dropouts we could 
have expected to hnd under the assumption that the quasar 
helds are “random” and show no enhancement due to the ex¬ 
istence of the quasars. The best comparison datasets which go 
deep enough over a wide area are the Hubble Space Telescope 
ACS imaging of the GOODS regions (Giavalisco et al. 2004) 
and the Subaru Deep Field (Kashikawa et al. 2004). These 
observations give a surface density of objects with z' < 25.5 
and i'-z' > 1.5 of 0.01-0.02 arcmin"^ (Dickinson et al 2004; 
Bouwens et al. 2004; Nagao et al. 2004). The total sky area 
we have surveyed with GMOS-North is 82 arcmin^. There¬ 
fore on the basis of the GOODS and SDF observations we 
would expect w 1 - 2 /'-band dropouts in our total area. Our 
hnding of one dropout is entirely consistent with the expecta¬ 
tions for a blank held. 

The results presented in this section show that these quasar 
helds do not exhibit an excess of luminous companion galax¬ 
ies. The magnitude limit of z' = 25.5 corresponds to a UV 
luminosity of ±1500 = 2.5 x 10^®ergs“^Hz“* at a redshift of 
z = 6.3. This is equivalent to 2L* in the z « 6 galaxy lu¬ 
minosity function (Bunker et al. 2004) and an unobscured 
star formation rate of SFR = 3OM0yr“\ assuming the con¬ 


version given in Madau, Pozzetti <& Dickinson (1998). The 
few known galaxies at redshifts z ~ 6.6 discovered in narrow- 
band surveys have star formation rates derived from their UV 
luminosities comparable to this limit (Hu et al. 2002; Kodaira 
et al. 2003). For comparison, the millimeter detections of 
dust in SDSS J1048H-4637 and SDSS J1148H-5251 imply total 
star formation rates > lOOOM 0 yr“^ in the host galaxies of the 
quasars (Bertoldi et al. 2003a). 

A similar study to ours has recently been reported by 
Stiavelli et al. (2005). They have observed the quasar 
SDSS J1030H-0524 with ACS on the Hubble Space Tele¬ 
scope in the 1775 and zgso hlters. Their data reaches slightly 
deeper than ours and they identify seven objects with i-ns - 
Z 850 > 1.5 in the quasar held. Only 3% of randomly-selected 
areas in the GOODS helds show such a high density of 
red objects, indicating an excess in the held of this quasar. 
The single i'-z' >1.5 object we discovered in the held of 
SDSS J1030H-0524 which is plotted on Fig.l^is also found by 
Stiavelli et al. and a high-resolution image of it is shown in 
the upper panel of their Fig. 2. Four of the seven objects found 
with ACS are detected at ijj^ and have 1775 -zsso < 2. These 
objects are not red enough to be located at a redshift simi¬ 
lar to the quasar (Dickinson et al. 2004). A spectroscopic 
redshift has been obtained for one object and it reveals a red¬ 
shift of z = 5.970. The large redshift difference between the 
quasar and this galaxy means that if they are part of the same 
large-scale overdensity then the comoving scale of that struc¬ 
ture is larger than anything found in the low-redshift universe 
(Pandey & Bharadwaj 2005). 

6. CONSTRAINTS ON THE DARK MATTER HALO MASSES 

Early-type galaxies in rich clusters in the local universe 
show homogeneous, old stellar populations suggesting co¬ 
eval star formation at high redshifts (e.g. Lopez-Cruz, Bark- 
house & Yee 2004). Therefore, if the quasars (two of which 
have far-IR luminosities indicating vigorous star formation) 
are in proto-cluster environments one would expect to see 
massive companion galaxies forming stars. It has been sug¬ 
gested that there could be a suppression of star formation in 
the vicinity of a strong UV photo-ionizing held such as that 
of a quasar (Couchman & Rees 1986). This is because the 
radiation would heat the diffuse gas in the IGM to a tem¬ 
perature of T ~ lOOOK and inhibit the cooling process by 
which gas condenses to form stars. However, these effects are 
likely to be signihcant only in low mass halos and will not 
be such an issue for the more massive halos hosting galaxies 
with SFR > 3 OM 0 yr“' which we are sensitive to. 

The goal of our observations was to constrain the masses of 
the dark matter halos hosting these quasars by considering the 
clustering of star forming galaxies around them. In fact, we 
have been unable to identify any companion galaxies suggest¬ 
ing that these quasars may not be residing in the most massive 
halos at this epoch, in line with the dynamical mass measure¬ 
ment of the host galaxy of SDSS J1148H-5251 by Walter et al. 
(2004). However, given the large uncertainties in transform¬ 
ing from observable quantities such as UV star formation rate 
to typical Lyman break halo mass and duty cycle at z > 6, it 
is difficult to directly turn our non-detection into a robust up¬ 
per limit for the halo masses. In this Section we show that 
scatter in the Mbh—Mdm relation, combined with a steeply 
falling mass function, biases the host halo masses of the ob¬ 
served high-redshift quasars to lower values, suggesting that 
the inferred low host halo masses do not necessarily imply 
a break-down of the low redshift Mbh-Mdm correlation at 
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z > 6. 

An upper limit on the halo mass hosting luminous quasars 
can be calculated from a comparison of their space densities. 
This method assumes that quasars reside in the most massive 
halos which exist in sufficient quantity to host all the observed 
quasars. The nine luminous quasars at redshift 5.7 < z < 6.5 
discovered by the SDSS (which all have black holes with mass 
= 1 - 4 X 10® Mq via the Eddington argument) give a 
space density of (6.4 ±2.4) x 10“^®Mpc“^ (Fan et al. 2004)^. 
The halo mass function at z = 6.3 was determined from the 
mass function fit to the numerical simulations of Sheth & 
Tormen (1999). The quasar space density corresponds to the 
space density of all halos with Mdm > 1.5 X IO'^Mq. This 
limit is cosmology dependent, and is obviously strongly de¬ 
pendent on the high mass tail of the Sheth & Tormen mass 
function, which is only relatively sparsely sampled in the sim¬ 
ulations used to determine this fit. However, in support of this 
functional form, Barkana & Loeb (2004) show that the tail 
provides a good approximation to the numerical mass func¬ 
tion at extremely high redshifts, provided that the lack of large 
scale modes in the simulations is taken into account. 

An estimate of the host halo masses can be found by us¬ 
ing the estimated black hole masses and further assuming that 
there is a correlation between black hole and halo masses. 
Such a correlation has been shown to exist at low redshifts 
(Ferrarese 2002) via a combination of the well known stellar 
bulge-black hole mass correlation and a correlation between 
halo circular velocity and bulge velocity dispersion. The cor¬ 
relation determined by Ferrarese (2002) with a correction ap¬ 
plied by Bromley et al. (2004) is: 


10*^Mq 


~ 0.015 


IO'^Mq 


1.82 


( 1 ) 


The universality of this correlation is untested and in partic¬ 
ular whether it evolves with redshift is highly uncertain. We 
will refer to this non-evolving correlation between Mbh and 
Mdm as case A. 

Wyithe & Loeb (2003) have noted that in feedback- 
regulated black hole growth models (e.g. Silk & Rees 1998; 
Fabian 1999) the important parameter for the strength of the 
dark matter halo potential is the circular velocity, not the 
mass. The relationship between circular velocity and mass 
is redshift-dependent due to the universal expansion: Vc oc 
Mq'^( 1-l-z)*^^ (e.g Barkana & Loeb 2001). The feedback- 
regulated model of Wyithe & Loeb® leads to a redshift- 
dependent correlation between Mbh Mdm of 


Mbh 


io^Mq 

Since the dormant black holes observed in the local universe 
typically formed at redshifts 1 < z < 2 (e.g. Yu & Tremaine 
2002 ) we fix the normalizing factor / so that the correlation 
roughly matches that of Bromley et al. (2004) at z = 1.5. This 
condition is met by / = 0.0018. We will refer to this evolving 
correlation between Mbh and Mdm as case B. 


/ Mdm 


5/3 


(1+z)'/^ 


( 2 ) 


^ We neglect the likely existence of quasars just as powerful as those in the 
SDSS which have their UV luminosities substantially decreased by dust. The 
space density of obscured quasars at lower redshifts is comparable to that of 
unobscured quasars (Zheng et al. 2004). However the ratio of obscured to 
unobscured quasars at z = 6 is so uncertain that we do not attempt to make 
any correction to the quasar space density. 

* We omit the parameter [^(z)]^^® from the equation given in Wyithe & 
Loeb (2003) since it is close to unity and negligible compared with the un¬ 
certainties involved. 


Since the two possible scenarios are normalized to each 
other at z = 1.5, at higher redshifts the host halo for a black 
hole of a given mass will be less massive for case B than case 
A. Specifically, for a quasar at z = 6.3 with Mbh = 2 x 10® Mq, 
then case A gives Mdm = 5.2 x 10 ^^Mq and case B gives 
Mdm = 1.3 X IO^^Mq. A simple comparison of the case A es¬ 
timate with the space-density limit shows that there are insuf¬ 
ficient halos with masses of Mdm = 5 x IO^^Mq at these high 
redshifts to host the number of observed quasars. The sim¬ 
ple conclusion to draw is that the correlation between black 
hole mass and halo mass does not remain the same at z > 6 
as it is locally. A similar conclusion was reached by Brom¬ 
ley et al. (2004) by considering the total black hole mass 
density built up by quasars. Wyithe & Loeb (2005) showed 
that the evolution in the 2dF quasar correlation function over 
the (admittedly narrow) range 1 < z < 2 was more consistent 
with the evolving Mbh “ Mdm correlation than the redshift- 
independent one. In what follows we will use both case A 
and B correlations, whilst remembering that prior evidence 
supports case B over A. 

An important ingredient in the application of such correla¬ 
tions that is often overlooked is the effect of intrinsic scatter 
in the relation. We now incorporate the scatter in these cor¬ 
relations to estimate the expected host halo masses. We will 
show that including the scatter makes a significant difference 
due to the steepness of the halo mass function for these rare 
halos. This is because lower mass halos are much more abun¬ 
dant, so low mass halos which are outliers in the Mbh-Mdm 
correlation could be more abundant than high mass halos on 
the correlation. 

We begin with the dark matter halo mass function fit of 
Sheth & Tormen (1999). This is evaluated at a redshift of 
z = 6.3. Above a mass of ^ IO^^Mq, the mass function is ex¬ 
ponentially declining. The black hole mass within each halo 
is given by the case A or B correlations. We now add scatter 
to these black hole masses. The intrinsic scatter in the corre¬ 
lation between black hole mass and stellar velocity dispersion 
at low redshift is in the range 0.25-0.30 dex (Tremaine et 
al. 2002). The conversion from stellar velocity dispersion to 
halo mass is likely to add in considerably more scatter, so we 
consider 0.4 dex to be the minimum amount of scatter in this 
correlation at low redshift. This scatter may increase substan¬ 
tially at higher redshifts, particularly for galaxies which are in 
the process of growing their black holes or accreting matter 
from the surrounding IGM. To bracket the uncertainty in the 
scatter we consider four different cases with 0.4, 0.6, 0.8 and 
1.0 dex scatter in the correlation. Although the shape of the 
scatter is unknown even locally, we will adopt a lognormal 
distribution as the simplest shape. 

By identifying the most luminous z ~ 6 quasars across a 
large fraction of the sky, the SDSS has selected the most 
massive, active black holes at this epoch (under the assump¬ 
tion that luminous quasars at such early times are accret¬ 
ing at the Eddington limit). Therefore we take a black hole 
mass cut of Mbh > 10® Mq and then calculate the probabil¬ 
ity distribution of dark matter halos which host such massive 
black holes given the correlations A and B and the various 
amounts of scatter in the correlation. The results are shown 
in Fig.|3 Including scatter has considerably decreased the ex¬ 
pected masses of the halos hosting these quasars. Even for 
the case with minimum scatter there is a shift of ~ 0.5 in 
logjgMDM- For the maximum scatter case considered, this 
shift in logjgMDM is ~ 1.5. 

Walter et al. (2004) use resolved CO emission to measure 
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log M„„ / Mo log M„„ / Ma 


Fig. 5.— Probability distribution for the mass of the dark matter halo host¬ 
ing a black hole with mass >10® Mq at z = 6.3. These curves were generated 
assuming correlations between Mbh and 3 Tdm of the form A and B (see text) 
with lognormal intrinsic dispersion cr. The vertical dotted (dashed) line shows 
the case A (B) Mdm expected for Mbh = 2 x 10® Mq without considering the 
combined effects of scatter and the steepness of the halo mass function. These 
plots clearly show that the expected halo masses are considerably lower than 
obtained by simply using the Mbh “Mdm correlation without scatter. 


a dynamical mass for the inner 2.5 kpc of SDSS J1148-1-5251 
to be ~ 5 X IO'^Mq. They argued this to be at least a factor 
of 10 lower than expected by naive extrapolation of the low 
redshift correlation between black hole mass and stellar ve¬ 
locity dispersion. We have shown here that using the case B 
redshift-dependent correlation and/or including the effects of 
a reasonable amount of scatter (tT> 0.6) lead to an expected 
halo mass at least an order of magnitude lower than using the 
naive correlation of case A. This provides a simple explana¬ 
tion for the observations of Walter et al. and those presented 
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in this paper, without the need for the Mbh-ZWom correlation 
to break down at high redshift. A consequence of this is that 
the SDSS quasars may not be identifying the most overdense 
regions of the high-z universe as is commonly assumed. 

7. SUMMARY 

We have presented deep imaging in the i' and z' biters of 
the belds of three of the most distant known quasars. For 
two of the three quasars, there are foreground galaxies at 
projected distances of a few arcsec. However, the faintness 
of these galaxies suggests they have relatively low masses 
which, combined with the large impact parameters, provide 
only a weak gravitational lensing magnibcation that would 
not signibcantly alter the derived quasar luminosities. 

A search for star forming galaxies at the redshifts of the 
quasars was carried out. No enhancement in Lyman-break 
galaxies above that in random belds is found. The UV con¬ 
tinuum star formation limit reached is 3OM0yr“', compara¬ 
ble with known z w 6.5 galaxies. The lack of companion 
galaxies could be indicative of these quasars residing in lower 
mass dark matter halos than previously believed. We consider 
the effect of scatter in correlations between black hole mass 
and halo mass and bnd that this scatter means that we expect 
the most massive black holes to be found in a range of halo 
masses, typically much lower than found by use of a simple 
correlation between the two quantities. 

To obtain a more accurate estimate of any overdensity 
of star forming companions will require much deeper data. 
Whilst few-orbit ACS observations are able to bnd a few pos¬ 
sible companions (Stiavelli et al. 2005), a proper census re¬ 
quires something of similar depth to the Hubble Ultra Deep 
Field (reaching 0.1 L* in the z = 6 galaxy luminosity func¬ 
tion with a limiting magnitude of z' = 29; Bunker et al. 2004) 
and such observations will therefore require a large increase 
in sensitivity such as from the James Webb Space Telescope. 

Thanks to the Gemini-North queue observers for execut¬ 
ing the observations presented here. We thank Stuart Wyithe 
and Meghan Gray for interesting discussions. Thanks to the 
anonymous referee for some useful suggestions to improve 
the paper. 
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